The superconducting state in the B2H6 compound at 360 GPa 
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In the paper, the thermodynamic properties of the superconducting state in the B2H6 compound 
have been characterized. The pressure of 360 GPa has been taken into account. The calculations 
have been carried out in the framework of the Eliashberg formalism for the wide range of the Coulomb 
pseudopotential: ^* £ (0.1, 0.3). It has been found that the critical temperature (Tc) varies in the 
range from 147 K to 87 K, depending on the assumed value of the Coulomb pseudopotential. The 
ratio of the energy gap to the critical temperature (_Ra = 2A (0) /fcsTc) significantly exceeds the 
value predicted by the BCS theory: f?A £ (4.24,3.98). In the similar manner behaves the ratio 
of the specific heat jump to the heat of the normal state {Rc = AC (Tc) /C^ (Tc)), namely: 
Rc e (2.33, 2.17). The parameter R„ = TcC^ {Tc) / (0), where He (0) is the thermodynamic 
critical field, ranges from 0.144 to 0.168. 

PACS numbers: 74.20. Fg, 74.25.Bt, 74.62. Fj 
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The intensive researches on the properties of the su- 

STConducting state have lasted for over a hundred years 
. Their main goal is to get the material in which the 
superconducting state would exist at the room tempera- 
ture. So far, the highest critical temperature {Tq) equal 
to 164 K has been measured in the HgBa2Ca2Cu308+y 
superconductor, located under the pressure at ~ 31 GPa 
0. However, due to the lack of the acceptable theory 
of the superconducting state in cuprates, it is difficult to 
answer the question whether in the present group of the 
materials might be possible to get even higher value of 
Tc i- 

In 1968, Ashcroft has suggested that the metalhc hy- 
drogen under the action of the high pressure would be 
the room-temperature superconductor [J. Later, after 
the complicated numerical calculations, it has been found 
that the pressure of the metallization (pm) for the hydro- 
gen molecular phase equals about 400 GPa Thus, in 
the questioned pressure range, the existence of the super- 
conducting state with the very high critical temperature 
can be expected [1] , , [l] , • 

Increasing in the value of the pressure up to ^ 500 GPa 
causes the dissociation of the hydrogen into the atomic 
phase. It should be noted that near the pressure of the 
molecular dissociation {p = 539 GPa), the critical tem- 
perature can be equal to 360 K [Toj . 

Currently, the highest value of Tc for the metallic hy- 
drogen is expected for the pressure at 2 TPa. In the 
considered case, the critical temperature can reach up 

600 K [m, [13, [i3- 

Due to the fact that the superconducting state in the 



metallic hydrogen can be formed only for very high pres- 
sures, scientists began to look for other physical system, 
in which Tc takes high value, while Pm is relatively low. 
Based on the survey, it has been found that the most 
interesting group is the family of the hydrogen-rich com- 



pounds [IJ, [15[, [161, PJ: M, lli, |2Q|- For exam- 
ple, in the silicon compounds Si2H6 [p = 275 GPa) and 
SiH4 (H2)2 [p = 250 GPa), the maximum of the criti- 
cal temperature takes the values of 173 K and 130 K, 
respectively [2l|, [l^l- 

In the present paper, we have described the results 
obtained for the superconducting state in the B2H6 com- 
pound {p = 360 GPa) [11]. Due to the high value of the 
electron-phonon coupling constant (A = 1.32), the cal- 
culations have been carried out in the framework of the 
Eliashberg formalism It should be noted that this 
formalism is the natural generalization of the classical 
BCS theory [25 1. 

On the imaginary axis ( j = \/— T), the Eliashberg equa- 
tions create the infinite system of the non- linear algebraic 
equations with the integral kernel, which allows to deter- 
mine the order parameter function (0„ = (j)(iu)n)) and 
the wave function renormalization factor {Z^ = Z (iw„)). 
The quantity ujn denotes the n-th Matsubara frequency: 
LOn = (7r//3)(2n- 1), where [3 = (ksTy^ {ks is the 
Boltzmann constant). The order parameter is defined by 
the ratio: A„ = (/)„/Z„. 

The open form of the Eliashberg equations can be writ- 
ten as: 



TT ^ A(iw„ - iw,,,.) - /.i*6' (wc - |w„i|) 

^ m=-M 
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A {iLUn - ii^m) 



iZm: (2) 
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where the pairing kernel for the electron-phonon interac- 
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tion is given by the formula: 



A(z) 



dn- 



(3) 



The spectral function {a^F (£7)) for the B2H6 compound 
under the pressure at 360 GPa has been calculated in the 
paper [23j . The maximum phonon frequency is equal to 
368.5 meV. 

The parameter ^* is called the Coulomb pseudopoten- 
tial and serves for the modeling of the depairing Coulomb 
correlations [2^; 9 denotes the Heaviside unit function, 
and cUc is the cut-off frequency; Wc = 30„iax- 

When carrying out the numerical calculations the finite 
number of the equations should be taken into account. 
It turns out that above simplification does not make the 
significant error to the final result, if the value of the 
considered temperature is not too low. In the study, we 



have adopted M = 1100 and T e (Tq = 25K, Tc), which 
ensured the convergence of the functions and Z„. 

It is worth noting that the Eliashberg equations have 
been solved using the numerical procedures, which have 
been tested and discussed in the papers [13] and (28j . 

From the physical point of view, the Eliashberg equa- 
tions on the imaginary axis allow to determine the critical 
temperature, the free energy, the thermodynamic crit- 
ical field, and the specific heat of the superconducting 
state. However, with their help it is impossible to ac- 
curately calculate the physical value of the energy gap 
(2 A) and the electron effective mass (m*). To estimate 
the parameters 2 A and m*, the solutions of the Eliash- 
berg equations should be analytically continued from the 
imaginary axis to the real axis (w). Such is the purpose 
of the Eliashberg equations in the mixed representation 
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M 
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where (a;) and / (cj) denote the Bose-Einstein and 
Fermi-Dirac functions, respectively. 

Note that the Eliashberg equations in the mixed rep- 
resentation have been solved using the pro cedures dis- 
cussed and tested in the papers [30| and [3l|. 

In the first step, basing on the numerical analysis, the 
possible range of the critical temperature has been de- 
termined. It has been found that Tq varies in the range 
from 147 K to 87 K for pi* € (0.1,0.3). Thus, regardless 
of the physical value of the parameter /i*, in the B2llg 



compound, the occurrence of the high temperature su- 
perconducting state is highly expected. 

It is worth noting that the value of the critical temper- 
ature can be determined by the McMillan or Allen-Dynes 
formula [H], [s^l- However, the obtained results are sig- 
nificantly depressed in comparison to the results obtained 
directly from the Eliashberg equations. 

The forms of the order parameter and the wave func- 
tion renormalization factor on the imaginary axis have 
been presented in Fig.[TJ Basing on the obtained re- 
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FIG. 1: (A)-(C) The order parameter and (D)-(F) the wave function renormahzation factor for the selected values of the 
temperature and the Coulomb pseudopotential. The first 200 values of Am and Zm have been presented. 



suits, it has been found that the order parameter values 
strongly decrease with the increasing temperature and 
the Coulomb pseudopotential. In turn, the wave func- 
tion renormahzation factor is much less dependent on T 
and /.J*. 




50 100 
T K 



150 





A □ 
A '-' 






A □ 






□ 

□ 




i: 

□ 


□ 

□ 

f 






/o 
/ o 

y o o 


n*=o.i 

H*=0.2 
H*=0.3 




A 


50 


100 


150 



TK 



FIG. 2: (A) The order parameter Am=i and (B) the renor- 
mahzation factor Zm=i as a function of the temperature for 
selected values of the Coulomb pseudopotential. The sym- 
bols represent the exact numerical results. The solid lines 
have been plotted on the basis of the formulas (6) and (7). 

The complete dependence of the order parameter 
and the renormahzation factor on the temperature and 
the Coulomb pseudopotential can be most conveniently 
traced after plotting the functions /S.rn=i {T , fx*) and 
Zm=i {T,fi*). The obtained results have been shown in 
Fig. [2 

It should be noted that from the physical point of 



view, the quantity 2Am=i (T, /i*) approximately deter- 
mines the value of the energy gap at the Fermi level. 
On the other hand, the function Zm=i {T,ij,*) is related 
to the effective mass of the electron. In particular: 
TO* ~ Zm=i {T, /J,*) nie, where the symbol TOg represents 
the electron band mass. 

Let us note that the numerical results presented in 
Fig. [2] can be parametrized by the means of the following 
formulas: 



T 



and 

Zrn^l(T,fJ.*) = Zrn=l{^J.*) 



(6) 



(7) 



T 



where F = 3.75. The functions A^^i (/z*) and Zm=i {n*) 
can be written as: 



A^^i (^*) = 108.6 {y.*Y - lOO.V + 35.8, (8) 



and 



Zrn=i (a^*) = -0.484 (/i*)^ + 0.467Ai* + 2.163. (9) 

The value [Tc] need to be calculated on the basis 

of the expression: Zm=i {Tc) = 1 + A = 2.32, where 



On the basis of the solutions of the Eliashberg equa- 
tions, the free energy difference between the supercon- 
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ducting and normal state has been calculated [34 



AF 



27r 
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E 

ri=l 



(10) 
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The symbol p (0) denotes the value of the electron density 
of states at the Fermi level; and are the wave 
function renormalization factors for the superconducting 
(S) and the normal state {N), respectively. 

The obtained results have been shown in the lower 
panels in Fig.[3l It is easy to note that the in- 
crease of the Coulomb pseudopotential results in the 
strong decrease of the free energy value. In particular: 
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FIG. 3: (Lower panels) The dependence of the free energy 
on the temperature for the selected values of the Coulomb 
pseudopotential. (Upper panels) The thermodynamic critical 
field as a function of the temperature. 

In the next step, the thermodynamic critical field has 
been determined: 



He 



V-Stt [AF/piO)]. 



(11) 



The results have been shown in the upper panels in Fig. 
[3l The destructive impact of the Coulomb pseudopoten- 
tial on the value of the thermodynamic critical field has 
been determined as: [He (0)]^*^o 3 / i-Hc (0)] 
0.55, where i/c (0) = -ffc (To). 

The specific heat of the superconducting state (C^) 
has been calculated on the basis of the formula: 



J^*=o.i 



C = C^ + AC. 



(12) 



where the specific heat for the normal state (C^) is the 

e^ fx") 

linear function of the temperature: ^Bpio) ^ ^' ^^'^ 
Sommerfeld constant is equal to: 7 



kBpiO) 

; |7r2(l + A). The 



difference in the specific heat between the superconduct- 
ing and normal state has been determined using the for- 
mula: 



AC (T) 
kBp{0) 



lf[AF/pm 

P d{kBTf 



(13) 
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FIG. 4: The specific heat of the superconducting state and 
the normal state as a function of the temperature for the 
selected values of the Coulomb pseudopotential. 

Fig. m presents the dependence of specific heats on the 
temperature and the Coulomb pseudopotential. It can 
be seen that at the critical temperature exists a char- 
acteristic jump, which is marked by the vertical dashed 
line. In addition, it should be also noted that the value 
of the specific heat jump strongly decreases with the in- 
creasing of Coulomb pseudopotential. In particular, the 
ratio [AC(rc)]^.^„.3/ [AC(rc)]^*=o.i equals 0.55. 

Basing on the obtained results, the values of the two 
characteristic dimensionless ratios have been calculated: 



Rh 



c) 



(0) 



and 



Rc 



AC(rc) 

(Tc) ■ 



(14) 



It has been found that in the range of the considered 
values of the Coulomb pseudopotential, the parameter 
Rh increases {Rh S (0.144,0.168)) with the increas- 
ing value of p.*, whereas the parameter Rc decreases 
{Rc G (2.33,2.17)). It should be boldly underlined that 
in the framework of the BCS theory, the ratios Rh and 



Rc arc the universal constants: [Rh 



BCS 



0.168 and 



[^cIbcs ~ 1-4^ The difference between the Eliash- 
berg predictions and the BCS theory results from the 
existence of the strong-coupling and retardation effects, 
which appear in the B2H6 compound. 

By using the functions (/)„ and one can solve the 
Eliashberg equations in the mixed representation. The 
obtained results for the order parameter have been shown 
in Fig.[SJ 
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FIG. 5: The order parameter on the real axis for the selected values of the temperature and the Coulomb pseudopotential. 
Additionally, the rescaled Eliashberg function has been plotted {70a'^ F (fl)) . 
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FIG. 6: The order parameter on the complex plane for the selected values of the temperature and the Coulomb pseudopotential. 
The lines with the symbols have been obtained for lj £ (0, fimax); the lines without the symbols correspond to the range of the 
frequencies from fimax to cjc. 



It has been found that for the low frequencies the non- 
zero value takes only the real part of the function A (a;). 
This result indicates no damping effects [IHl • Let us note 
that the imaginary part of the order parameter becomes 
non-zero only for the higher frequencies. 

In addition, it can be easily seen that the functions 
Re [A (w)] and Im [A (uj)] have a rather smooth course. 
However, the existing weak maximums or minimums are 
related to the corresponding peaks in the Eliashberg func- 
tion. 

It also seems worthy to plot the values of the order pa- 
rameter on the complex plane (see Fig.[n|). In the present 
case, the characteristic shapes of the deformed spirals can 



be observed. It should be noted that the radius of the 
considered curves becomes strongly reduced with the in- 
creasing temperature and the Coulomb pseudopotential. 

The curves presented in Fig. [5] allow in an easy way 
to characterize the effective potential for the electron- 
electron interaction. In the case when the curve meets the 
condition Re [w] > 0, the potential is pairing Thus, 
on the basis of Fig. [HI it has been found that the effective 
interaction between electrons leads to the formation of 
the superconducting condensate in the frequencies range 
from to ^ 0.90riinax, if ^* =0.1. The increase in the 
value of the Coulomb pseudopotential causes the nar- 
rowing of the considered range of lj from the side of the 
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higher frequencies. 

The energy gap at the Fermi level has been determined 
on the basis of the equation: 



A(T) = Re[A (w = A(r))], 



(15) 



From the physical point of view, the most interesting 
is the value for the lowest temperature: 2A(0), where 
A (0) = A (To). As a result of the calculations the fol- 
lowing has been obtained: 2A (0) £ (26.77, 14.91) meV 
for ^* e (0.1,0.3). 
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FIG. 7: (A) The wave function renormalization factor on the 
real axis. Additionally, the rescaled Eliashberg function has 
been plotted (3a^F(r2)). (B) The wave function renormaliza- 
tion factor on the complex plane. The lines with the symbols 
have been obtained for lj £ (0, ilmax); the lines without the 
symbols correspond to the range of the frequencies from f^max 

to LJc- 



The wave function renormalization factor on the real 
axis for T = Tc and n* — 0.1 has been shown in Fig. 
[7] (A). Next, we have presented the courses of the renor- 
malization factor on the complex plane for the selected 
values of the temperature. 

The second solution of the Eliashberg equations al- 
lowed us to calculate the exact value of the electron effec- 
tive mass: m* = Z (w = 0)me. As a result, it has been 
found that in the entire temperature range, from Tq to Tc 
the electron effective mass is large and it reaches its max- 



imum value at the critical temperature {[r, 



2.56). 



In summary, in the paper we have determined all rele- 
vant thermodynamic parameters of the superconducting 
state in the B2H6 compound. The pressure of 360 GPa 
has been taken into account. 

It has been found that, regardless assumed value of 
the Coulomb pseudopotential, the critical temperature is 
high. In particular: Tc G (147,87) K, for /i* changing in 
the range from 0.1 to 0.3. 

Other thermodynamic parameters differ significantly 
from the predictions of the classical BCS theory. This 
is especially visible for the low values of the Coulomb 
pseudopotential. In particular, it has been proven that 
the dimensionless ratios, characterizing the derogation 
from the results of the BCS theory, take the values: Ra € 
(4.24,3.98), Rc € (2.33,2.17), and Rh e (0.144,0.168). 



On the basis of the results above, it is very easy to 
estimate the value of the dimensionless parameter: 



Ra 



2A(0) 
ksTc ' 



(16) 



The obtained result has the form: Ra e (4.24,3.98). 
It should be noted that in the framework of the BCS 
theory the ratio Ra is the universal constant and it is 
equal to 3.53 [2^. Hence, it can be easily noticed that 
the BCS model is too simple to predict in the correct way 
the physical value of Ra for the B2H6 compound under 
the pressure at 360 GPa. 
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